Chiral networks in twisted graphene bilayers under interlayer bias 
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A twisted graphene bilayer exhibits a triangular Moire pattern in the local stacking, that smoothly 
alternates between the three basic types AA', AB' and BA'. Under an interlayer bias U, the latter 
two types develop a spectral gap, characterised by opposite valley Chern numbers. We show that 
for large enough Moire periods and bias (angles smaller than ^0.2° at U ^ 90 meV) these regions 
become depleted electronically, and topologically protected chiral modes appear at their boundaries. 
This gives rise to a delocalised chiral network of the Chalker-Coddington type, composed of valley 
current vortices. Simultaneously, a discrete set of localised states at resonant energies develop. 
Clear signatures of this exotic electronic state are predicted to arise in the spectrum and the optical 
conductivity. 



In spite of graphene being a zero gap semiconductor [T], 
the practical difficulty to induce a finite gap has hindered 
its potential for nanoelectronics applications. In this re- 
gard, graphene bilayers were shown early on to possess 
one significant practical advantage over their monolayer 
counterpart. For Bernal-type stacking of the two lay- 
ers (also known as AB'or BA' stacking [2 ), a band gap 
can be opened around each valley by applying an electric 
bias between the two layers [3l H] . Gaps up to 250 meV 
have been achieved experimentally OE]. This property 
allows a graphene bilayer to be electrically tuned into a 
very high mobility semiconductor [7], greatly extending 
its potential for applications. 

The gap tunability of graphene bilayers is intriguing 
also from a fundamental point of view. Firstly, it depends 
critically on the precise AB' or BA' interlayer stacking. 
An AA' bilayer, for example, is a good metal for any rea- 
sonable interlayer bias. Secondly, the gap can be made 
to "change sign" by inverting the bias U into —U. Per- 
forming such an inversion between two adjacent regions 
in space gives rise to two topologically protected chiral 
(TPC) modes per valley, that flow without resistance 
along the interface between the two regions [8 [9 . Anal- 
ogous TPC channels arise in a range of systems, in which 
the gapped local band structure changes topology across 
a boundary, e.g. at the edges of Hall bars, or the surface 
of topological insulators [TOl [11] . Lastly, and rather sur- 
prisingly, a change of the topology of a tuned gap can be 
induced also with a uniform bias U by smoothly transi- 
tioning from AB' to BA' stackings [12 , which are related 
by mirror symmetry. This can be achieved by applying a 
different strain to each layer, and also gives rise to TPC 
modes. 

In this work we address the problem of the formation 
of TPC modes in graphene bilayers with a generic (non- 
minimal) stacking. These are known as twisted graphene 
bilayers, and are characterised by a finite interlayer rota- 
tion angle. We find that TPC modes may arise in these 
systems without strain under a uniform interlayer bias 



U. These chiral modes arrange into a triangular network 
of valley-polarized current running along zero mass lines 
[13] . This is a variation of what is commonly known as 
a Chalker-Coddington network p^HTH] , Superimposed 
onto this delocalized network, a discrete set of states 
strongly localized away from AB'/BA' regions emerge at 
resonant energies. Such restructuring of electronic states 
occurs for energies |e| < \U\/2, and produces clearly mea- 
surable signatures in the density of states (DOS) and the 
optical conductivity o-xx{oj). Interestingly, the mecha- 
nism that gives rise to the chiral network also appears, 
although for different reasons, in graphene monolayers on 
boron-nitride (h-BN) substrates [T31IT7HT9]. The lattice 
mismatch between graphene and h-BN, however, does 
not allow for an arbitrarily large Moire period, so that 
unlike in twisted bilayers, the chiral network is not ex- 
pected to fully develop in that case. 

Band topology and chiral modes. — Let us briefly re- 
view some concepts on the band topology of graphene 
bilayers. Any electronic band defined in a finite Brillouin 
zone (BZ) has an associated integer topological invari- 
ant (TTJ [2Ql[22] , known as Chern number C, which de- 
termines e.g. its contribution to the Hall conductivity 
o'xy = Ce'^/h when the band is completely filled [23j- 
27 . In 2+1 dimensions, the Chern number is defined 
in terms of eigenstates \ip{k)) as C = ^ J^^(PkF{k), 
where F{k) = 2lin{dk^ip{k)\dkyip{k)) is the Berry curva- 
ture. A biased AB '-stacked graphene bilayer develops a 
gap around its two valleys, located at points K and K' 
in the BZ. The Berry curvature of the valence band is 
peaked around these two points and has opposite sign, 
so that C = 0. However, a (non-quantized) valley Chern 
number may be defined by confining the k integral to the 
region around a single valley, C^o = ^ JkC) d'^kF{k). 
If the bias is small as compared to the interlayer cou- 
pling, U < j±, Ck = —Ck' ^ sign(/7) is approximately 
quantised. This valley Chern number plays then the role 
of a proper topological charge [28] , as long as valleys are 
kept decoupled by any perturbation (i.e. as long as the 



bilayer, and U{r) itself, are smooth on the scale of the lat- 
tice constant ao). In particular, it implies the emergence 
of two TPC modes per valley and spin along an interface 
where U (and hence Cj^o) smoothly changes sign [8]. 
Likewise, since a BA'-stacking has opposite valley Chern 
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numbers than AB' (C^f/^ 



-C^^) ), a smooth transition 



between these two stackings will also produce two TPC 
modes per valley and spin [12 . This result is a version 
of the bulk-boundary correspondence, which states that 
N = \C2—Ci I mod 2 TPC modes emerge where the Chern 
number changes from Ci to C2. The mod 2 is absent in 
our case {N = 2) since we consider only valley-preserving 
perturbations, hence "weak" topological protection, see 
Ref. [11] A transition between stackings arises naturally 
in a bilayer with a relative rotation between layers, the 
so-called twisted graphene bilayer. 

Continuum theory of twisted bilayers. — A twisted 
graphene bilayer is parameterised by the relative rota- 
tion angle 6>, with = denoting the AA' bilayer. Ge- 
ometrically, such rotation produces a periodic triangular 
Moire pattern in the local stacking, which smoothly al- 
ternates between the three minimal types AA', AB' and 
BA'. For an angle of rotation < 30°, the Moire pat- 
tern has a period Lm = ao/ [2sm{0/2)] (ao = 2.46A is 
the monolayer Bravais period). Electronically, the low 
energy sector of the system does not depend on the de- 
tailed crystallography [29H331 . and is well described by 
a continuum theory in which the two valleys are decou- 
pled, and the Dirac points in each layer are offset by 
a momentum AK = 2|i^| sin(6>/2), with \K\ = 47r/3ao. 
In the absence of interlayer coupling, the two cones in- 
tersect around the M point of the superlattice, at energy 
eM{0) = hvFAK/2 (which is ^ 88meV x 6>(deg) for smah 
angles). This is the natural energy scale in the low en- 
ergy sector. The coupled bilayer is described within each 
valley by the continuum spin-degenerate Hamiltonian 
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which is written in the A, 5, A\ B' basis. Here n± = kx^ 
i{ky T Ai^/2), and the function V{r) = ^7^ ^.^^^^ e'^i" 
describes the periodic spatial variation of the interlayer 
coupling, with 7^ ^ 0.33eV. The Moire lattice vec- 
tors are denoted by ai^2 (so Lm = 1 0^1,2!), and conju- 
gate momenta are ^1,2, such that Qi • aj = 27r6ij (we 
also define gs = 0). The AA' sublattice is centered 
at the origin, with the AB'/BA' sublattices offset by 
±Sr = (ai - a2)/3. 

The continuum model has two distinct regimes, de- 
pending on the value of dimensionless parameter a{6) = 
\i^/^m{0). The regime of large angles (6> ::^ 1°, a <C 1) 
is amenable to perturbation theory in a, and has been 
extensively studied. It is characterised by the forma- 
tion of a van-Hove singularity [29l I34H37] at energy 
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FIG. 1. (Color online) (a) Bandstructure and (b) density of 
states (DOS) for a ^ = 4° twisted bilayer with and without 
interlayer bias U . At such angles, a bias merely shifts the two 
Dirac cones, creating a flat DOS plateau for |e| < \U\/2 [to 
the left of the dashed line in (b)]. 



CvH ~ eM(l — 2a) close to the M point [see peak in the 
DOS at ~ O.Gcm in Fig. [lib)] and a suppression of the 
Fermi velocity Vp ~ '^f(1 — 9a^) [29l[35]. In the presence 
of a finite bias /7, the band structure in the large angle 
regime reacts in a predictable fashion, merely shifting 
the two Dirac cones by ±/7/2, see Fig. [ija), since their 
eigenstates mostly reside in different layers in this case. 
This results in a low energy density of states that is flat 
instead of linear for |e| < |/7|/2, see Fig. [ijb). The 
DOS profile is reminiscent of the perfect AA'-stacked bi- 
layer's, whose band structure is also composed of two 
shifted Dirac cones. Above |/7|/2, the DOS is largely 
unaffected by the bias. 

The regime of small angles has been less explored and 
has a much more complex structure. In the unbiased 
case, secondary van Hove singularities above the first ap- 
pear and move lower in energy as 6 is decreased. As 
each singularity approaches zero energy, it morphs into a 
quasi-flat miniband, which is accompanied by a vanish- 
ing Fermi velocity Vp. This occurs at specific ("magic") 
twist angles Oi^ the highest being Oi ^ 1.05° (which cor- 
responds to Lm ^ 52ao) |38||39j. In general however, 
this simple picture can only make sense of the spectrum 
for angles down to around ~ 0.35°. Below this 6>, 
the [7 = electronic structure becomes rather intricate. 
With the notable exception of a resilient AA'-localized 
state [40l |4T] pinned at around e ~ 0.1 7eM (which is 
evolved from the original large angle van-Hove singular- 
ity), the DOS profile, as shown in the bottom curve of 
Fig. |2|a), exhibits no clear structure (although see Ref. 
39), and states are typically delocalised (Fig. |2|f)). 

The small-angle electronic structure becomes much 
simpler as U is increased. The AB' and BA' regions de- 
velop a local gap, with opposite valley Chern number. If 
U is above certain threshold, and these regions are large 
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FIG. 2. (Color online) (a) Density of states of a Lm — 290ao {0 ~ 0.2°) twisted bilayer with increasing bias U . The region 
above the dotted line denotes the chiral network regime, (b-f) Periodic spatial profile of electron probability density for several 
characteristic eigenstates in the spectrum [arrows in (a)], all chosen at one of the K points. Black colour corresponds to zero 
density, and white to maximum. Pannel (b) shows a fully developed chiral network (Chalker-Coddington) state, with arrows 
indicating the orientation of the valley currents. Pannels (c), (d) and (e) show various resonant states, all localized away from 
AB'/BA' regions by the finite bias. Pannel (f) shows a [/ = delocalized state. 



enough (small enough angles), the bias tends to electron- 
ically deplete them. As a consequence, strong confine- 
ment of states away from the depleted AB'/BA' regions 
becomes possible. A discrete set of sharp localized reso- 
nances appear in the spectrum, as seen in the top curve 
of Fig. [2|a). These localized states arrange spatially in 
a variety of patterns [see Fig. [2|c,d,e)], but always away 
from the AB'/BA' regions. For off-resonant energies, the 
DOS exhibits a weak and featureless background that is 
in marked contrast to the intricate small- angle, [/ = 
DOS profile, and to the Dirac-like linear background of 
the large- angle DOS. This weak background corresponds 
to delocalized TPC modes running along AB'/BA' inter- 
faces, which develop due to the spatial modulation in the 
valley Chern number of the gapped regions. The collec- 
tion of chiral channels arrange in a triangular Chalker- 
Coddington chiral network, as shown in Fig. [2|b). Each 
channel carries opposite current for opposite valleys. The 
pattern of valley currents for a given sign of U is fixed by 
the band topology, and is shown by white arrows. 

This new electronic phase requires strong depletion in 
AB' and BA' regions to fully develop [13 , which in turn 
requires large bias and small angles. One can approx- 
imately quantify this condition by considering the sim- 
plified one-dimensional problem of an abrupt AB'/BA' 
interface with a uniform bias. Using the two component 
approximation (valid for [/ < 7^), an interface mode de- 
cays like - e-^l^l, with A = ^/TPyZ/2vF [8 • We look 
for the value of U that results in an amplitude decay 
of the interface mode of at least, say, 95%, within the 
radius of the Moire BA' region of the twisted bilayer. 



R = Lm/'^V^- This is half the distance between AB' 
and BA' region centers, see Fig. [2|b). Inserting x = R 
this results in the following condition for the development 
of the chiral network 
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where p = —4^/3 log 0.05 = 20.76, and l± = Hvf/j^ ^ 
7.3ao is the interlayer coupling length [42 . 

This analysis suggests that for a realistic value of 
/7 ~ 90 meV, the network requires a Moire period 
Lm > 290 ao {0 < 0.2°), the value chosen in Fig. [2| 
For this angle, the range of U and e that satisfies the 
criterion has been marked by a dotted boundary in Fig. 
|2|a). We see the boundary roughly correlates with the 
crossover into the chiral network regime, characterised by 
sharp DOS peaks on a weak smooth background. 

The DOS depends only on the density of energy eigen- 
values of the system. Other equilibrium observables, such 
as the optical conductivity axx{^)^ contain additional in- 
formation related to matrix elements between different 
eigestates. It can be computed using the Kubo formula 
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where F(e, e') = Tr [^(e - H)jJ{e' - H)j,i and j, is the 
current operator in the x direction. Quite famously, the 
low-frequency limit of the optical conductivity is univer- 
sal in a neutral graphene monolayer (Jxx{^ ^ 0) = ctq = 

2 
^ ^ , which corresponds to a universal infrared transmis- 
sion at normal incidence T = (1 + 7ra/2)~^ ^ 97.7%, in 
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FIG. 3. (Color online) (a) Zero-temperature optical conduc- 
tivity Re[cr(c(;)] at neutrality (/i = 0) for increasing bias U 
and twist angle like in Fig. [2] It exhibits complex fluctua- 
tions around 2ao = 7re^//i at small U (dashed line). As the 
system transitions into the chiral network phase (above thick 
dotted line), the fluctuations are transformed into a smooth, 
non-monotonous background corresponding to transitions be- 
tween chiral states, with weak superimposed resonances from 
transitions between localised states. Some of the latter [ar- 
rows in (a)] are shown in (b) for U — 180 meV and increasing 
frequency [arrows from bottom to top]. 



terms of the fine structure constant a [43l|44]. In an unbi- 
ased, large-angle twisted bilayer, axx{^) exhibits features 
associated to the van-Hove singularities J451447] . As the 
system enters the low-angle regime, the optical conduc- 
tivity exhibits a complexity similar to the DOS, in the 
form of fluctuations around a rough average of 2ao (which 
corresponds to two decoupled monolayers). Under an 
increasing interlayer bias, however, axx{^) becomes re- 
markably simple, like the DOS, see Fig. [3J It evolves into 
a smooth background that is non-monotonous in cj, and 
which falls below 2cro, representing the optical absorp- 
tion of the delocalised chiral network. This background 
exhibits superimposed peaks that reflect transitions be- 
tween sharp localised states [see arrows in panels (a) and 
(b)]. Such peaks, however, are much suppressed as com- 
pared to those in the DOS, due to the small overlap of 
localised eigenfunctions under the action of the current 
operator jx. Hence, absorption due to the chiral network 
becomes the dominant feature of (Jxx{^) at low frequen- 
cies. For frequencies larger than |/7|/2, complex features 
are still present, of a magnitude similar to those of the 
unbiased case. 

In conclusion, we have shown than an interlayer bias 
generates a starkly different low energy electronic struc- 
ture for Bernal-stacked, AA'-stacked and twisted bilayer 
graphene. While the former are semiconducting and 
metallic, respectively, the latter develops a novel spectral 



structure at low twist angles and realistic values of the 
bias, as a consequence of the depletion of AB' and BA'- 
type regions in the Moire pattern. The periodic modu- 
lation of the local band topology gives rise to the forma- 
tion of a delocalized chiral network of protected states of 
the Chalker-Coddington type, whose density of states is 
smooth and decreasing up to energy ±\U\/2. Localized 
states also appear at specific resonant energies in various 
spatial configurations away from the depleted AB'/BA' 
regions. These features are visible both in the spectral 
density and in the optical conductivity, which become 
strongly modified as the bias exceeds the threshold for 
the formation of the chiral network phase. 
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